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Abstract: This paper considers the construction of 3-band compactly supported symmetric biorthog- 
onal wavelets. Given an interpolatory refinable function, we present an explicit iterative 
algorithm for constructing its dual scaling function with desired regularity, and discuss the 
properties of the constructed scaling function, such as its symmetry and decay. Finally, 
several design examples are given. 
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1 = Introduction 


In [1], Daubechies constructed a family of compactly supported orthonormal scaling func- 
tions and their corresponding orthonormal wavelets with dilation factor 2. Daubechies has 
proved that there exists no dyadic compactly supported orthogonal scaling function which is 
symmetric, except for the Haar scaling function. For the 3-band case, it is known from [2,3] 
that there also exists no compactly supported symmetric orthogonal scaling function which 
is interpolatory. Hence, such scaling functions have to be biorthogonal. In [4], Ji and Shen 
constructed a class of compactly supported biorthogonal wavelets and orthogonal interpolatory 
wavelets with dilation factor 2. Methods for constructing M-band biorthogonal wavelets were 
also given in [5,6]. Examples of biorthogonal B-spline type wavelets with certain regularities 
were obtained in [7]. In [8], the authors considered the design of 3-channal adaptive biorthogo- 
nal filter banks via lifting. In this paper, we will consider the construction of 3-band compactly 
supported symmetric biorthogonal wavelets, i.e., given a interpolatory refinable function, we 
present an explicit iterative algorithm for constructing its dual scaling function with desired 
regularity. 

Recall that a multiresolution analysis (MRA) with dilation factor M , İs a Sequence of nested 
subspaces of L?(R), i.e., ++- C V1 C Vo C Vi C---, where 


Vj = span{f(MJz — k) : k € Z} 


for some f(x) € L?(R), and 


Uv=P(R), fv = {0}. 


jEZ jEZ 
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The function f(x) is called the scaling function. Let f(x) be a compactly supported scaling 
function with dilation factor M. Then f, r (x) #0, and f(z) satisfies the refinement equation 


f(a) = SJ af(Mz—k), Soa =M, (1) 
keZ keZ 
where cx are real, and ck # 0 for only finitely many k € Z. Define a trigonometric polynomial 


P(z) = >> ckz", ze, 


M fet 
2 Preliminary results 


Define masks as 
mo(w) = Vi au hret, olw) = D2 GeT., (2) 


Here, we assume that only finitely many hn, gn are nonzero. Some of our results can be extended 
to infinite sequences that have sufficient decay for |n| — oo. Then we define ¢, ¢ by 


= Il mo(M-iw), $w) = Il Mmo(M iw). (3) 


j=l j=l 


These infinite products can converge only if mo(0) = 79(0) = 1, i.e., if 
2O e a 
n n 


If (4) is satisfied, then the infinite products in (3) converge uniformly and absolutely on 
compact sets. So ¢(w) and (w) are well defined C% functions. What is more important, if 
the infinite products in (3) converge in L?(R), then two scaling functions ¢(x), ¢(x) will be 
well defined according to [9]. From Lemma 3.1 in [10], ¢(z) and (x) have compact support. 

Let ¢(x), (x) € L?(R) be the scaling functions, and Mo(w), Mmolw) be their masks, re- 
spectively. ġ(x) is called the dual function of olx) if 


(p(x), d(x = k)) = 6,05 kez. 


If (x£) and ¢(z) are biorthogonal scaling functions, then their masks must satisfy 


ail 


Y mofu + 2+ BE) = 6) 


Definition 2.1 A continuous scaling function ¢(x) is interpolatory if it satisfies ¢(k) = 
k,05 kez. 
From [4], it is easy to know that a compactly supported continuous function @(x) is inter- 
polatory if and only if 
So dw + 2kr) = 1. 


kez 
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Let mo(w) be its mask, then it satisfies 


AE 2rk 
X m(w+ ©) =1, Vwe R. (6) 
M 
k=0 
Condition (6) is called the interpolatory condition. Condition (6) is only a necessary condition. 
A necessary and sufficient condition in terms of the refinement mask is given in [11]. From [4], 
the autocorrelation of an orthonormal scaling function is interpolatory. 


3 Dual functions of interpolatory scaling functions 


In this section we only consider the case of M = 3. Here we consider 


(7) 


dee gay 
— -N 
P(w) = (E) Qu), 
where 
k . 
Q(z) = D a;27, a-k =a 
j=H-k 
and all a; are real. 
If the mask P(w) is of the form in (7), then it is said to have N-regularity. It is easy to 
know that the Fourier transform of ¢(z) is given by 


d(w) = [| Pw). (8) 
j=l 


Theorem 3.1 Let P(w) be a 27-periodic function of the form in (7) that satisfies the 
interpolatory condition (6). Define 


H(w) := P?(w) + 3P(w)(1— P(w)) + 2P(w+ *) P(w i =). (9) 


Then the function P(w)H(w) satisfies the interpolatory condition (6) and the coefficients of 
H(w) are all real. Also P(w)H(w) is of the form in (7). 
Proof First let H(w) be the function defined by (9). Then 


P(w)H(w) = 3P?(w) — 2P3(w) + 2P(w)P(w $ p(w + =), 


3 3 
P(w+ =) (w+ =) = 3P?(w+ =) — 2P3(w + =) 
+2P(w)P(w+ =)\P(w+ =), 
P(w+ 2) a (w+ =) = 3P?(w+ =) — 2P3(w + =) 


+2P(w)P(w + 2) P(w $ =). 
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Therefore, 
S P(w + a) y(u + aá = |Pw) +P(w + =) + P(w + a =]; 
k=0 


Next, from (7) and (9), to prove the coefficients of H(w) are real, it suffices to show the 
coefficients of 


out Belu 3) 
are real. Observe that 


Since 


eter ete) =2 o 1, Vj} nEzZ. 


Hence the coefficients of (10) are real. 

Since P(w) is symmetric about the origin, consequently, P(w)H(w) is of the form in (7). 
So Theorem 3.1 holds. 

Let ¢(x) be the symmetric interpolatory scaling function whose mask P(w) is defined by 
(7). We can use Theorem 3.1 to construct the dual function of ¢(z). First define the dual mask 
by 


P(w) = P2(w) + BP(w)(b ~ Plw)) + 2P (w+ =) P(w n = . 


Lemma 3.2 Suppose that the function ¢(x) defined by 


$w) = [J #3 


j=l 


is in L?(R). Then ġ(x) is a dual function of d(x). 
Proof The proof is similar to the Theorem 2.2 in [4]. 
Theorem 3.3 Let the mask P(w) be of the form in (7), and Q(z) satisfy 


2k 
Qloz) #0, Q(0)=1, w= T, k=l 
If there exists some n = 0 such that 
B, = max |Q(w)Q(w/3)-- Q(w/3")| < alae (11) 


then (w) satisfies the decay condition 


à -1/2- i 1 logB 
< 1/2-€ Z aA n 
\o(w)| < CU+ \w]) , with e=N 5S ees 


> 0, (12) 


where C > 0 is a constant. What is more, $(x) € L?(R) and satisfies the refinement equation. 
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Proof From the Appendix B in [6], it is easy to know that if the condition (10) can be 
satisfied, then (11) holds consequently. So ¢(x) belongs to L?(R). By applying the Theorem 
2.17 in [12], there exists a function ¢ € L?(R) and it satisfies the following refinement equation 


PO PE Ener) = deo, 
j=l 


Corollary 3.4 Let P(w) be a mask of the form in (7) and aj, j = 0,1,--- ,k be all 
nonnegative. If the coefficients a; satisfy 


k 
5 aj < 3N 1/2, 
j=0 
then there exists a scaling function ¢(z) € L? (R) corresponding to P(w) defined by 
$(w) = J | P(w/3*). 
j=l 


Now, we give an iterative method. 
Iterative construction Let ¢(x) be the given interpolatory refinable function whose 
mask P(w) is given by (7). Let Po = P(w). For k =1,2,--- , then 


(i) define 
Pp := Phot [PR (w) + 3Px-1(w)(1 — Pk-1(w)) + 2Pp—1 (w + T) Pe (w + EN]: 
(ii) define 
P, = a 


(iii) define 


From the above theorem and algorithm, if in the k-th step, the conditions in Theorem 3.3 
are satisfied, then we can iterate these steps to get a function o(w) with better regularity. 


4 Examples 


In what following, we give two examples of biorthogonal 3-band refinable functions. 
Let Po be the mask of the scaling function $(x), which is symmetric and interpolatory, 
where 


Po(w) = yo(Ltetey? 


First, let k = 1. By using the iterative construction, we get 


~ L42422\2/—427! 411 — 42 
PT eine a ue nese es pt ee ae 
PLY ( 3 ) ( 3 ). 
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It is easy to verify that P, satisfies the conditions in Theorem 3.3. So it can generate a dual 
function ¢1(x) of ¢(x). Clearly, PoP, is the mask of another interpolatory scaling function 
p(x). Similarly, we can get a dual function (xz) of y(x). 

Let k = 2. By using the above algorithm, we get P», which also satisfies the condition in 
Theorem 3.3. So it can generate a dual function ¢2(x) of ¢(x). Let 


~ 1 
P = 3 3 Gnz”. 
We obtain that 


[gn]! 17 = [-2.9725e-006, —5.9446c-006, — 8.9153e-006, 5.5141e-005, 0.00035997, 
0.00060456, — 0.00069907, — 0.0065288, — 0.013667, 0.015167, 0.053461, 
0.079964, — 0.097909, — 0.25377, — 0.28033, 0.32117, 0.9687, 1.4269, 
0.9687, 0.32117, — 0.28033, — 0.25377, — 0.097909, 0.079964, 0.053461, 
0.015167, — 0.013667, — 0.0065288, — 0.00069907, 0.00060456, 
0.00035997, 5.5141e-005, — 8.9153e-006, — 5.9446c-006, — 2.9725e-006]. 


The following figures show the graphs of (x) and the dual scaling functions of ¢(#). From the 
graphs we know that the larger of the k, the smoother of the dual function. 


os s| 
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Figure 1: Graph of $(x) Figure 2: Graph of ¢1(2) 


G 
2 — 12 —_— g r a ai 2 3 r -— 
' 


a a a a o a a ‘e a M Aa s a as o s a s a as OO a a ‘a ‘o a a s 


Figure 3: Graph of $2(z) Figure 4: Graph of (x) Figure 5: Graph of (x) 
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